Spin Density Wave Fluctuations and p-wave Pairing in Sr2Ru04 
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Recently a debate has arisen over which of the two distinct parts of the Fermi surface of Sr2Ru04, is the active 
part for the chiral p-wave superconductivity. Early theories proposed p-wave pairing on the two dimensional 
7-band, while a recent proposal focuses on the one dimensional (a, P) bands whose nesting pockets are the 
source of the strong incommensurate spin density wave (SDW) fluctuations. We apply a renormalization group 
theory to study quasi-one dimensional repulsive Hubbard chains and explain the form of SDW fluctuations, 
reconciling the absence of long range order with their nesting Fermi surface. The mutual exclusion of p-wave 
pairing and SDW fluctuations in repulsive Hubbard chains favors the assignment of the two dimensional 7-band 
as the source of p-wave pairing. 
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While there is a general consensus in favor of a chiral p- 
wave superconducting (SC) state in Sr2Ru04 [1-4] stabilized 
predominantly in either the 2-dimensional (2D) 7 band, or the 
pair of quasi-lD (a, /3) bands, a debate has arisen recently 
over which of the two is the active source of the supercon- 
ductivity. The debate has been triggered by the failure to ob- 
serve the persistent edge currents associated with the chirality, 
though this also has been controversial [3, 5, 6]. Microscopic 
derivations of a chiral p-wave pairing state in the 2D-scenario 
have been proposed based on a 2D Hubbard model by a T- 
matrix approach [7], third perturbation theory [8] and func- 
tional renormalization group (RG) calculations [9]. On the 
other hand, Raghu et. a/. [5, 10] have recently argued for the 
quasi- ID-scenario as more compatible with the missing edge 
currents and have provided a microscopic justification for it by 
using a RG theory, but only in the limit of weak interactions. 

A closely related and competing phenomenon in Sr2Ru04 
is the strong spin density wave (SDW) fluctuations at an in- 
commensurate nesting wavevector spanning the Fermi sur- 
faces of the 4/3 filled (a,/3) bands [11]. SDW fluctuations 
at this wavevector Q = (27r/3, 2tt/3) [12] were recently re- 
ported at room temperature and at energies as high as 80 meV 
[13]. The SDW peaks at Q which combine nesting in both 
nearly ID Fermi surfaces, grow as the temperature (T) is re- 
duced and saturate at the crossover to 3D Fermi liquid behav- 
ior at T-sD 60 K [1 1] when the resistivity starts to show a 
behavior [14]. Smaller peaks were observed at the wavevec- 
tors (7r,27r/3) and (27r/3,7r). To date these have been dis- 
cussed within a random phase approximation (RPA) scheme 
[15-17, 30]. Because of the highly nesting character of the 
(a, /?) Fermi surface, it is necessary to choose a very weak in- 
teraction in RPA with a value typically an order of magnitude 
smaller than standard estimates [18, 19]. 

In this Letter we show that treating the ID character of the 
(a, P) bands in a RG scheme, can explain the strong SDW 
fluctuation and reconcile the absence of the SDW long range 
order at T > T^o using a standard value for the interac- 



tion as observed in Sr2Ru04. Furthermore, our RG scheme 
shows mutual exclusion of p-wave pairing and SDW fluctu- 
ations in repulsive Hubbard chains, and a sharp suppression 
of the SDW fluctuations at low frequency in the p-wave SC 
state. Such a suppression has not been observed in previous 
neutron experiments and will be a challenge to explain within 
the quasi-lD-scenario. 

We start from a ID RG [20] treatment of the single chains 
in the (a, /3) bands. This includes the important cancellation 
between particle-hole and particle-particle graphs, which is 
absent in RPA. The properties of single chains in one-loop 
RG were derived in an early application of RG to condensed 
matter systems [20]. With repulsive interactions, the SDW 
and triplet superconductivity (TS) response functions have 
a power law form with divergences to infinity or zero, as 
T ^ Q. The phases with enhanced SDW and suppressed 
TS and vice versa, are separated by a Quantum Critical Point 
where the exponent 6 changes sign. In the one-loop approx- 
imation, 6* = 92 — .91/2, where gi (32) is the dimensionless 
scattering coefficient for backward (forward) processes. The 
response functions in ID RG behave differently from that in 
RPA, which gives a finite scale divergence in the SDW re- 
sponse. This RPA scale corresponds to a large temperature 
(up to lO'' K) in Sr2Ru04 if we use typical values of the intra- 
chain hopping t = 0.3 eV [2] and Hubbard repulsion U — 2.2 
eV [18, 19]. 

To begin with, we consider the hybridization and interac- 
tion between the two orbitals, dxz and dyz, which give rise to 
perpendicular chains. To obtain estimates of the mean field 
transition temperature to SDW order in Sr2Ru04 we use the 
results of numerical calculations [21, 22] to extend the one- 
loop RG calculations to stronger interactions. Next, we ex- 
amine the interactions between parallel chains. This leads to 
a theory, which gives a strong SDW fluctuation without long 
range order at low temperatures. 

The effect of the hybridization and spin-orbit coupling 
among the two orbitals can be described by the following 
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Hamiltonian [5] 
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Here ct is the electronic creation operator in the orbital 

/c,m,cr 

TO with momentum fc and spin cr, and the angular momentum 
operators and spin operators are represented in terms of the 
totally anti-symmetric tensor = itamn and Pauli matri- 
ces a, respectively. For the system we are interested in, the 
strengths of the mixing and spin-orbit coupling are t" w O.lt 
and r/ w O.li, respectively [23, 24]. 

When the above perturbation 811 is taken into account, the 
quasi-particle spectrum opens a gap 2A near (zb/ci?, ±fci?) with 
A = ■\/ (3i"/2)2 + if-. Therefore, the dispersion for the dxz 
orbital is modified to 

eg ~ WF(|fcx| - kp) + sgn(|fc:j| - kp)L{\ky\)\ (2) 

with L{x) the Lorentzian function centered at kp. The 
Green's function describing quasi-particle excitations in the 
dxz orbital is Gxzikxi ky, icj„) = l/(icj„ — e^), and a similar 
counterpart in the dyz orbital. Then the bare SDW response 
function reads 



/Cxz ( 
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with N the total number of sites, and we work in units where 
the lattice spacing is unity. For q = {2kp, qy), we have, after 
a analytic continuation to real frequency with w = 0, 
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with Eq = At the band width and vp the Fermi velocity, where 
we have set ft = 1. The bare SDW can be approximated by 
the following interpolation scheme, in the Umit of small A, 



X™(q,T) = x™(q,T)«^ln^. 

nvp rjQ 



(5) 



A standard RG calculation, by including the particle- 
particle and particle-hole graphs, gives the dressed suscepti- 
bility in terms of intraorbital interactions: 
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with = U /2'KVp in the Hubbard model. To deal with the 
case of strong interactions, we compared our results with the 



T-dependence of spin susceptibility from the Monte Carlo 
method[21, 22], and found that 9 is screened to be Q.W, which 
we denote as 6* from now on. From the expression above, one 
can see that due to the hopping between the two orbitals, a fi- 
nite low energy cut-off A appears, killing the divergence as 
0. 

Next we introduce the inter-orbital interactions between 
the two 4(i orbitals, with the generic form of the on-site 
interaction[18] 



Hi = Y^ Un 



[U'nrnanna + {U'-JH)nrnanna 



ni<n 

(7) 

with U' the interorbital Coulomb repulsion and the Hund's 
rule coupling, in addition to the intraorbital repulsion U. Here 
cl^g. is the creation operator for electrons in the orbital m with 
spin cr at the same site, and Uma is the corresponding number 
operator Here we note that the intraorbital interaction U has 
been taken into account in the exponent of Eq. 6. 

To incorporate the above interactions, we define the joint 
response function for SDW, by including the orbital indices 



XH{q,i^)^- e'''^TrOiq,T)OHq,0))dT, (8) 
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where 



i„J-yU,,ntir)], (9) 



with q = {2kp,qy) and {qx,2kp) for dxz and dyz orbitals, 
respectively. 

To first order in the inter-orbital interactions, the only non- 
vanishing term is 
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Gyz{k' + Q, IL0'„ + in)Gyzik',lUj'J, (10) 



with a corresponding diagram shown in Fig. 1 (a). Note that 
the wave vector for the response function Q — {2kp, 2kp) is 
the same for both dxz and dyz orbitals due to the conserva- 
tion of momentum in the scattering process in Fig. 1. Another 
important consequence is that only the Hund's rule coupling 
contributes to the SDW response function, while other onsite 
interaction terms in Eq. 7 are not involved. This result origi- 
nates from the spin configuration in Fig. 1. In this sense, the 
Hund's rule coupling assists the spin-flip processes between 
different orbitals. An intuitive physical picture is that the spin- 
flip processes are coherent even in different orbitals, due the 
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FIG. 1 . (Color online) (a) The lowest order Feynman diagram for the 
spin-spin correlation function connecting propagators from different 
orbitals via the Hund's rule coupling. The red and blue lines stand 
for electrons in the xz and yz orbitals, respectively. The solid and 
dashed lines correspond to electrons belonging to the branches con- 
taining +kp and —kp in the one-dimensional model, respectively. 
Due to the conservation of momentum at the vertex, Q is locked to 
be {2kp, Ikp). (b) SDW vs. paramagnetic phase diagram in the pa- 
rameter space of A and Jh at T = from Eq. 12. The estimated pa- 
rameter region for Sr2Ru04 is indicated in the paramagnetic phase. 
In this region, Jh ranges from 0.13 eV to 0.4 eV, the renormalized 
and bare values of the Hund's rule coupling. 

ferromagnetic Hund's rule coupling between the two orbitals. 
Dynamical mean-field theory found that the Hund's rule cou- 
pling is important in Sr2Ru04 [18]. In our calculations below 
we use a Gutzwiller renormalized values of Jh ~ 0.13 eV 
to take into account the strong onsite repulsion between holes 
[25]. 

The full dressed joint SDW response function in Eq. 8 is 
obtained by first including the intra-orbital interaction U in 
a RG scheme, which means that the bare bubbles in Fig. 1 
(a) are replaced with the dressed ones in Eq. 6. Due to the 
absence of the low energy divergence in Eq. 6, the Hund's rule 
coupling Jh can be treated in a RPA-like method, leading to 

The divergence of xh{Q, T) in Eq. 11 gives estimate of the 
mean field transition temperature to long range SDW order. 



The first term gives an upper limit on T^^w Hund's 
rule coupling, which is about 50 K, similar to T^d- The pres- 
ence of the second term, of order 10"^ K, guarantees that the 
ground state is paramagnetic. To illustrate the underlying 
physics, we construct the phase diagram in the Jh ~ ^ pa- 
rameter space as shown in Fig. 1 (b). The parameter set for 
Sr2Ru04 is in the paramagnetic region, but near the phase 
boundary with SDW fluctuations. Therefore, the strongly en- 
hanced response function by the Hund's rule coupling at Q 
naturally explains the strong enhancement of the SDW signal 
near Q in the experiments [11]. 

Next, we briefly study how the interchain tunnelling be- 
tween parallel chains affects the SDW in the quasi-one di- 
mensional (a, (3) bands of Sr2Ru04. The interchain hopping 
t± alone would give rise to a singular SDW response function 
at the wave vector (2fci?,7r) since the tight-binding approxima- 
tion preserves the perfect nesting for quasi-one dimensional 
systems [20, 26]. However, in the case of Sr2Ru04, the Fermi 
surfaces the (a, (3) bands are distorted due to hybridization 
and spin-orbit coupling between orbitals, as discussed previ- 
ously. Therefore, the nesting property at {2kp,TT) is lost, and a 
strong enhancement for the SDW fluctuation is not expected. 

Another mechanism to affect the SDW response function 
at {2kp,TT) is the superexchange interaction Jex = 4t']_/U be- 
tween two neighboring parallel chains. But a rough estimation 
yields Jex ^ Jh, since t± is only about 0.026 eV [2]. Com- 
bining the above two effects for the parallel chains, the spin 
fluctuation response at {2kF, tt) should be much weaker than 
that at {2kF, 2kF), as observed in the experiment [13]. 

Finally, we consider the effects due to the SC pairing in 
ID bands on the magnetic response. Earlier measurements 
did not show a change in the magnetic response at Q upon 
cooHng through the SC transition at T^{= 1.5 K) [1 1]. 

To make our analysis more transparent, we restrict the dis- 
cussion to one dimension, and thus consider the following 
Hamiltonian H = Hq + Hint, where Hi„i — U rii^riii 
is the Hubbard onsite interaction term, which can be reduced 
to the standard form describing different scattering processes 
with .91 = (72 = U/ttvf [20, 27]. To incorporate the SC pairing, 
we consider Hq as 

-ff = ^ VF{\k\ - kF)c\^Cka 

+ ^[A(fc)4/_,^ + H.c.], (13) 

k 

which models a one-dimensional electrons with p-wave SC 
pairing A(fc) and can be solved exactly in the mean-field ap- 
proximation. The assumption that spin-orbit coupling locks 
the d vector along the crystal c axis has been made, consistent 
with the polarized-neutron scattering experiment in Sr2Ru04 
[28]. 

A powerful tool to study the properties of _ffo is the ap- 
plication of the normal and anomalous Green's functions 
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Go-cr' {k, iojn) and F^^i {k, ia;„) [29], which read 



and 
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Here = VF{\k\ — kp) and Ao-o-'(fc) = A{k)azi<Jy. Near 
the Fermi surface, we have, due to the odd parity, A(fc) = 
sgn(fc) Aq, with Aq the SC gap near the Fermi surface. 

Interestingly, after some algebra, we find that the particle - 
particle and particle-hole diagrams just differ by a minus sign, 
similar to the case without SC pairing. For example, the 
particle-hole bubble diagram can be expressed as 



T 
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To the leading order in the logarithmic accuracy, this expres- 
sion is reduced to 
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with Q(x) the Heaviside function, and we have performed an 
analytic continuation to real frequency uj at zero temperature. 
The structure of this expression is also similar to its counter- 
part in the normal state. And due to this analogy, the RG flow 
equations for the interaction constants gi and 52 should be the 
same as the non-SC case [20, 27]. Therefore, in the case of 
6* > 0, the SDW phase is expected as usual. 

A standard RG analysis yields the final results for the re- 
sponse function in the SC state as follows: 



Rex'^(2fcF,c^) = 
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As is shown in Fig. 2 (b), if one looks at the low-energy 
properties cj 2Ao, the response function diverge as x 
\ui — 2Ao|^^ This result indicates that the transition to 
superconductivity in the ID bands will open a gap in the low 
energy spectra at wavevector Q [16, 30]. While early neutron 
scattering experiments by Braden and coworkers [11] did not 
show a change in low energy spectra at Q, a more complete 
investigation would be worthwhile to definitively decide if an 
SC gap opens up in the ID (a, (3) bands at the onset of super- 
conductivity at = 1.5 K. 

In summary, we have applied a RG scheme starting from 
the ID analysis for single chains, to explain the strong SDW 



G{-k,i.Sl-iui„) 
-S 5-> >■ 



ys/\^s/s/r ^- 




G(k.iu„) 

(a) 



G(fe+2*;jr,!ui„+iri) ; G(fc,ttj„) 



F(k,iu„) \ IF'' (k+2kp,iLj„+in) 



(c) 



(d) 



3 
2.5 

2 

1.5 









Imx^"^ 




, (e) 



0.5 1 1.5 2 2.5 3 3.5 4 
(iVAn 



FIG. 2. (Color online) (Upper panels) Structure of the vertex dia- 
grams, (a) and (b) are for the particle-particle channel, whereas (c) 
and (d) are for the particle-hole channel. The solid and dashed lines 
conespond to electrons belonging to the branches containing -\-kF 
and — /cf in the one-dimensional model, respectively. The wavy lines 
stand for bare onsite interactions. (Lower panel) Response function 
as a function of at T = is shown in (e). Here x^'"(2fcF,<^) is 
scaled by Xo'^ 



Rex (2fcF,tJ = 0), and 9* is about 0.41 in the 



system of interest. 



fluctuations and the absence of SDW order at temperature 
above the crossover to 3D Fermi liquid behavior with the 
strong onsite Hubbard repulsion estimated for Sr2Ru04. The 
mutual exclusion in ID RG theory of enhancement in the 
SDW and simultaneously in the p-wave pairing channel is in 
favor of the 2D 7-band as the source of the superconductivity. 
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